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1 Introduction 

o : 

. Matrix models appear again and again at the front line of theoretical physics, and every new generation of 
04 ' scientists discovers something new in this seemingly simple subject. It is getting more and more obvious that 
. matrix models capture the very essence of general quantum field theory and provide the crucial representative 
D ' example for the string theory [1]: matrix models play for string theory the role, which harmonic oscillator 
plays for quantum mechanics. Essential difference is, however, that we have a nearly exhaustive understanding 
CO ' of harmonic oscillator, while nothing like complete description of matrix models is yet available. Probably, the 
time is coming to begin a systematic analysis with the goal of building up a theory of matrix model partition 
' functions as the first special functions of string theory [2] . Various applications of matrix model techniques should 
. use these functions as building blocks for formulation of their results, thus separating the physical content of 
different applications from common mathematical formalism. The goal can be to build up an analogue of the 
. powerful free-field formalism, developed in 1980's, which allowed to reduce many problems in perturbative 
' string theory [3] (and related issues in representation theory of Kac-Moody algebras, 2d conformal field theory 
. and finite-zone solutions in integrable-systems theory) to almost classical set of special functions: Riemann's 
theta-functions [4, 5], associated with Riemann surfaces. The problem is that matrix model partition functions 
are non-trivial generalization of Riemann's theta-functions, which were not fully investigated in XlX-th century, 
therefore the theory should involve essential new ideas, both in physics and in mathematics. This is the reason 
_ , why the progress in this field - actually known under the nick-name of non-perturbative string theory 
- is considerably slower. More and more people begin to realize that generic problems of non-perturbative 
(~| , string theory, i.e. the general theory of phase transitions, have a lot in common, and these most interesting 
common properties are clearly represented at the simplest level of matrix models. Complications, introduced by 
sophisticated context of particular applications can and should be separated from important content, captured 
at the matrix model level. 

In development of matrix model theory one can distinguish between several stages. In fact, these stages are 
typical for every chapter of the string theory, and they characterize not so much the history, but the kind of 
questions that are posed and the level of abstraction in analysis of these questions. Whenever a new problem 
arises (or is re-addressed again), it is quite useful to understand its place in this hierarchy. 

2 Classical period: introduction of models 

During the classical period the main task was to study concrete phenomena, involving matrix models. Since 
in these notes we are not going to discuss applications, we just mention a few crucial theoretical methods, 
developed at this stage. 

The most important idea was to reduce the A^^-fold matrix integrals like 



X 



to TV-fold integrals over eigenvalues (pi of matrix-valued fields $ = C/^diag((/)i)C/. This idea, closely related in 
general context to the problem of gauge invariance, is technically based on the possibility to explicitly integrate 
over "angular variables" in the simplest situations, of which (1) is a representative example. The single known 
exactly solvable generalization, which goes slightly beyond this example, involves the Itzykson-Zuber integral 

[dU]exp{TTAU^BU), (2) 
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which is used in construction of two important classes of matrix theories: generahzed Kontsevich model [8, 9], 



Zgkm{L\W) - 
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(3) 



mxm 



and even further generalized Kazakov-Migdal-Kontsevich model [10, 11], 



ZG(KMy 



■ exp [tr (W„(M„) + M^ul^M^Uaf} + L^Ul^M^Uc,^'^ 



(4) 



defined for any graph with Hcrmitian fields Ma and background fields standing in the vertices and unitary 
fields Uai3 and Ua/j. on the links. In more general situations (sometime quite important for applications, e.g. 
involving non- trivial multi-brane backgrounds), angular integration is highly non- trivial and remains an open 
problem. The simplest of such difficulties are echoed in the Gribov copies problem in Yang-Mills theory, 
while treatment of more sophisticated unitary-matrix integrals face problems, typical for adequate treatment of 
quantum gravity. The class of theories where gauge (angular) variables can be integrated out in effective way 
is referred to as eigenvalue models [6] . Most of further development concerns eigenvalue models, which until 
recently remained implicitly synonymous to matrix models. 

The next idea of the classical period used the fact, that angular integration usually provides non-trivial, but 
very special measures on the space of eigenvalues, which is often made from Van-der- Monde determinants [12]. 
For example, (1) can be transformed into 



— 4>i)- Occurrence of determinants suggested two alternative technical methods of investigation of 
eigenvalue models: technique of orthogonal polynomials [13] and free-fermion representation [13, 14]. Occurrence 
of discriminants was not fully exploited yet, though it opens a set of interesting possibilities, both technically 
and conceptually. 

The last idea of the classical period which needs to be mentioned, is the idea of continuum limit, when 

the matrix size N tends to infinity. In fact there arc infinitely many different continuum limits relevant 
for different particular applications, - and only very few have been analyzed so far. However, from the very 
beginning the main fact was broadly realized: continuum limits of matrix models posses description in terms 
of Riemann surfaces [15, 16]. Today we know that occurence of spectral surfaces is a general phenomenon, and 
they occur already at finite values of [2]. 

3 First stringy period: generalizations and hunt for structures 

The first stringy period is characterized by the change of interests: from tools to theory. Instead of studying 
various applications and developing adequate technical methods to answer the problems, posed by these applica- 
tions, attention gets concentrated on the search and understanding of internal structures. Instead of "study" the 
main slogan becomes "deform and generalize" - this is the standard string theory method of revealing hidden 
structures. It is at this stage that the three main inter-related structures were discovered behind eigenvalue 
matrix models: rich Ward identities, integrability and CFT representations. 

3.1 Ward identities 

Occurence of Ward identities is the pertinent property of every integral (and thus of quantum mechanics 
and all its generalizations, like quantum field and string theory): they refiect invariance of the integral under 
the change of integration variable ~ an archetypical example of auxiliary field. However, this obvious hidden 
symmetry manifests itself in a rather sophisticated manner: as relation between various correlators in the theory. 
String theory normally deals with partition functions: the generating functions of all the correlators, summed 
up with the coefficients like tk in (1), which have the meaning of extra coupling constants, and can be considered 
as providing the deformation of the bare action. This formalism allows to treat Ward identities as equations 
for partition function, because the change of integration variables can be compensated by the change of the 
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coupling constants, if there are many enough [6]. In particular case of the integral (1) the equations are known 
as Virasoro constraints [17] or loop equations [15], 

L_{z)Z{t\N)=0, (6) 
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For the integral (3) the Ward identities can be written in two forms: either as a Gross-Newman equation [18], 

{W {d/dLtr) + L) Zgkm{L\W) = 0, (7) 

which generalizes a similar equation [6] for the Itzykson-Zuber integral (2), or as a set of peculiar W„-constraints 
(where n = dcg W'{z)) [19]. A particular case of these equations in the case of n = 2 is continuum Virasoro 
constraint [20, 21], 

= 0, m > -1, 
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where t/. =ri. + ^trL~'°/" and = J^^^ res (VF'(z))^^''^" dz = — |(5fc_3. See [6] for detailed discussion of Ward 
identities for eigenvalue models and [22] for more technicalities, related to generalized Kontsevich model. 



3.2 Integrability and RG evolution 

Integrability of matrix models moans that partition functions satisfy not only linear equations (Ward identi- 
ties), but also bilinear Hirota-type equations, 

A(Z(8)Z) = 0. 

Technically, the proofs rely upon determinant representations of eigenvalue models, which are, in turn, imme- 
diate corollaries of determinant structures of integration measures, see [6] and references therein. However, the 
true meaning of integrability remains obscure. Several inter-related ideas should be somehow unified to clarify 
the issue. 

First, integrability should express the fact that the system of Ward identities is rich enough: enough to 
specify the partition function almost unambiguously - up to some easily controllable degrees of freedom, a sort 
of zero modes with some clear cohomological interpretation. See [23] for the relevant notions of strong and week 
completeness. 

Second, bilinear equations normally involve the matrix size N in non-trivial way: for example, the simplest 
bilinear equation for partition function (1) - the lowest term of the Toda chain hierarchy - states [24] that 

This means that bilinear equations involve not only variations of the coupling constant tk, but also those of 
conjugate variables, like A^. However, in (1) only N - a. conjugate variable for to (in the sense that d log / dto = 
A'^) - is present, while there is nothing like conjugate parameters for all other i's. This can suggest that Hirota 
equations for such restricted partition function are too special to reveal their general structure. In particular, 
(9) holds only for the simplest phase of the model (1), other phases belong to multi-component generalizations 
of Toda chain hierarchy. 

Third, while coupling constants parameterize the bare action, i.e. the weight of summation over paths 

in functional integral, the conjugate variables should rather characterize the integration domain (range of 
integration, boundary conditions, the target space - whatever formalism one prefers to use in introducing 
them). In this sense conjugate variables are intimately related to generalized renormalization group flows, 
which are supposed to do exactly the same: describe the dependence of functional integral on the integration 
domain. 
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Relation between integrability and renormalization group is one of the main open problems in modern 
theoretical physics [1, 6, 25]. To get fully related, both these concepts should be considerably modified. To 

become a pertinent feature of all partition functions, integrability should not be restricted to ordinary t- 
functions [26], associated with free fermions and single- loop Kac- Moody algebras at level k = 1. Indeed, 
T-functions and bilinear Hirota-like equations can be defined for arbitrary Kac-Moody and even more general 
Lie algebras and quantum groups [27]. However, the theory of such r-functions is not reducible to that of 
Plucker (free-fermion) determinants and requires the full-scale application of the free-field formalism [3] and 
more sophisticated determinant formulas. Expressions of this type are expected to arise in the study of unitary 
matrix integrals, but only first steps are being done in this direction [28, 29]. 

Renormalization group theory also requires considerable generalizations, of which we put especially empha- 
size on two. First, as already mentioned, it should allow arbitrary variations of integration domain, not just a 
one-parameter cut-off procedure: renormalization group should study the changes of the shape of integration 
domain, not just of its volume. This seems to be already a widely accepted generalization. Second, renor- 
malization group theory should be made applicable to the study of fractal structures, which often reproduce 
themselves (exactly or approximately) at different scales. This means that the theory should allow non-trivial, 
periodic and perhaps even chaotic, renormalization group flows [30, 31], and this should not contradict the 
obvious uni-directional nature of these flows, generated by integrating out degrees of freedom. This entropy-like 
feature of renormalization group is usually expressed through the c-theorem [32], but in the case of systems 
with infinitely many degrees of freedom the c-function can easily have angular nature, and monotonic decrease 
of c does not contradict the existence of periodic motions. Besides the obvious example of self-similar fractals, 
today we already have some field-theory examples of such behaviour [33] - but only at the level of flows. It is 
not so easy to provide examples of partition functions (not just fiows) with such properties, and the reason for 
this is simple: we usually consider as solvable the systems where answers can be expressed through conventional 
special functions, and these - almost by definition - have oversimplified branching structure. We do not posses 
any language to describe more general - and more realistic situations. At the same time such language should 
exist, because what seems to be sophisticated (fractal and/or chaotic) phase structure is in fact a highly ordered 
and pure algebraic pattern [34] , not so much different from the algebraic-geometry backgroimd under the theory 
of conventional special functions. The jump between order (integrability) and chaos (infinite-tree structure of 
bifurcations or phase transitions) can be not so big, and most probably the careful study of partition functions 
already at the matrix-model level will reveal the deep inter-relation (duality) between the two concepts. 

This duality should play a big role in landscape theory, which studies the distributions of various algebro- 
geometric quantities on moduli spaces (either of coupling constants or on the dual space of different branches) and 
their interplay with renormalization group flows. The crucial hidden double-loop structure, which is supposed 
to become the basis of such theory [35], should be seen already at the level of matrix models. For old and new 
attempts to apply - the yet undeveloped - landscape theory to phenomenology see [1] and [36] respectively. 

3.3 Relation to conformal theory 

The third structure behind eigenvalue matrix models is that of 2d conformal field theories: partition functions 
of matrix models coincide with certain correlators in conformal models. This fact is of course related with 
the free-fermion and KP r-function representations of simplest partition functions, but it was quickly realized 
to have a more general meaning [37]: CFT representations are also closely related to Ward identities. In the 
simplest example of (6) the natural observation is that ^{d(j))'^{z) looks like the stress tensor T{z) of the 2d free 
scalar field, and therefore one can look for an intertwiner, relating operator L{z) in (6) with T{z) = ^{dip)'^{z) 
in the free field theory. The intertwiner in this case is just 

where (p{z) is the 2d scalar, (0| the left vacuum, annihilated by <f+{z), {0\(p+{z) = 0, and v{z) = J2k^kz''. 
Obviously, 

i(^)(o|ei'^(^)^'^(^) = (o|ei'''(^)^'^(^)i(a^)2(^), 

and acting on the right vacuum {dip)'^{z)\0) = 0, so that 

i_(z)(0|pi''(^)^^(^)Q|0) = 0. (10) 

Here Q is inserted in order to make matrix element non-trivial, and the only requirement is that Q commutes 
with T{z) = ^{d(p)'^{z). Such Q can be made out of the screening charges Q± = <^ e^^'^^^\ and Z{t\N) in (1) 
is reproduced if Q = Q^. Many generalizations are straightforward, among them the important class of quiver 
matrix models [37, 38] . An immediate desire is to find a representation of this kind for Kontsevich r-function 
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- solution of continuum Virasoro constraint (8), but this problem appears surprisingly complicated and still 
remains open, see [39, 40] for important but still insufficient steps towards its solution. 

4 Transcendental period: absolutization of structures 

The transcendental or second stringy period is characterized by absolutization of structures, revealed at the 
previous stage. This means that the logic is inverted: now the structures arc given, and wc look for generic 
objects, possessing such structures - expecting in advance that some of these objects can appear different from 
original matrix integrals. Today we are mostly at this stage, at least in the theory of matrix models, and only 
the first careful steps arc being done in attempts to understand the hidden meaning of the structures that we 
observed in applications - and continue to observe again and again in new physical systems when they come to 
our attention. 

Discussion of generic CFT correlators as well as of generic r-functions is old and important but still not 
a very constructive problem. It lies at the very core of landscape theory, but existing theoretical methods are 
too week to address the problem directly. The structure which can be rather effectively analyzed by available 
tools or by their straightforward modifications is that of Ward identities. Working on this structure one can 
also hope to develop stronger methods, applicable to analysis of the other two structures. To mention just one 
open question, illustrating the degree of ignorance in this field, it remains unclear to what extent the Virasoro 
constraints (6) and (8) per se, without explicit integral representations like (1), imply integrable structure of 
partition function. 

4.1 Ward identities 

Defining partition function as solution to Ward identities is a typical D-module-style definition, i.e. basically 
a problem from linear algebra, since Ward identities are linear equations (this is what makes examination of 
this structure much simpler than integrable one, related to quadratic equations). Since we do not expect any 
mysteries in linear algebra (or, better to say, we think that we already know all of them), the theory building 
is straightforward, but by no means trivial: we know what to search for, but exhaustive and explicit solution 
can be quite tedious. We are not yet at the stage when solution of a linear-algebra problem, especially infinite- 
dimensional, and all its properties can be predicted a priori. 

Normally, the set of questions one asks about the solution consists of several groups. 

4.1.1 Domain of definition 

First of all, one should decide what class of functions one wants the solution to belong to. Basically, there are 

two principal alternatives: formal series and globally defined functions. However, even if one chooses formal 
series, it is still necessary to specify, where singularities are allowed to occur. The simplest alternative is 
between series in positive powers of a variable (singularities are allowed at infinity of the Riemann sphere) and 
in all integer powers (singularities allowed at infinity and at zero). However, there arc many more possibilities: 
singularities can be allowed at other points (then negative powers of [z — a] can be allowed), the underlying 
bare Riemann surface need not be a sphere (then certain fractional powers of [z — a] are allowed) etc. More 
than this, different variables can have different kinds of singularities and we have infinitely many such variables. 
In different applications different requirements are imposed, and one and the same quantity - unique partition 
function = generic solution to Ward identities - can show up in absolutely different way: different breinches 
are relevant for different applications. What can look like a singularity from the point of view of particular 
application, can be nothing but a branching point or phase transition to another branch, which acquires 
a natural interpretation in terms of an absolutely different application. Once again this emphasizes that the 
problem of partition fimctions - even one particular partition function, one particular matrix model _D-module,- 
is almost undistinguishable from the central problem of the string theory [1], which can be formulated as a search 
of unique universal partition function (a universal object of quantum field theory), of which all other thinkable 
partition functions are particular branches and sub-families. Returning closer to the Earth, even speaking about 
formal series one - explicitly or implicitly - imposes certain requirements on allowed singularities of solutions, 
which often do not follow from the Ward identities themselves. Actually, the "rich" Ward identities, that we 
spoke about in the previous section, fix everything but the behaviour at singularities, they are rich enough to 
fully constrain the dependence on the choice of the action (on coupling constants), but not rich enough to fix 
conjugate dependencies, e.g. boundary conditions at possible singularities: in certain sense they fix exactly one 
half of possible freedom. This is a typical thing for the equation on a wave function in quantum mechanics to 
do (constrain gr-dependence but ignore p, or in other words, there are many solutions of the same Shroedinger 
equation, differing, say, by the values of energy): this observation is a starting point on the road, leading to 
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identification of Ward identities with Shroedinger-type equation and the partition function with a (string-field 
theory) wave function. 

Convenient way to speak about emerging freedom is in terms of "globally defined" functions: to fix it one 
can ask to solve Ward identities for functions on Riemann sphere or on torus or anywhere else. The problem 
is that in most interesting cases - and matrix model partition functions are not an exception - singularities 
are unavoidable. Moreover, small relaxations of compactness, like "functions with simple poles", which were 
enough for the formalism of free fields (and thus in dealing with arbitrary 2d conformal models), now are not 
sufficient: r-functions are different from conformal blocks because they have essential singularities, even if we 
take those defined on punctured Riemann surfaces of finite genus. Actually, things are even worse: interesting 
partition functions require this genus to be infinite, if at all defined. Therefore this language - though widely 
used - by itself does not provide unambiguous classification of solutions. 

4.1.2 Generating functions 

The next question - after the set of allowed functions is somehow specified - is how to represent the answer. 
There are few chances that the full answer will be any known special function: as already stated, the most 
important thing is that matrix models produce new special functions. What is important, however, these new 
functions are not too new, they are just one step forward more involved than the known special functions, 
namely the Riemann's theta- functions, which occupy the previous level of complexity. We expect that all kinds 
of simplifications of matrix model partition fimctions reduce them to quantities, already expressible through 
Riemann's theta functions. What can be these simplifications? Instead of considering a generating function 
of all correlators, one can select some one-parametric family, and such reduced generating function can be a 
candidate. One can take one or another limit (say, one or another continuum limit: naive, double-scaling, fixed 
genus, ...). One can take a ratio of different branches (a kind of monodromy matrix) or some more sophisticated 
combination of those. 

4.1.3 Different realizations 

As usual for a D-module, one can look for solutions of linear equations in different forms. 

Integral formulas like (1) are particular examples of one possible - integral - representation (or realization, 
to avoid confusion with group-theory representations) of solutions. For integral realizations the Ward identities 
play a role of Picard-Fucks equations: solutions are periods of the form which is converted into a full derivative 
by application of the corresponding operator (which generates Ward identity). Of course, in such situation one 
can integrate along any closed contour and different non-homological contours give rise to different solutions. 
Therefore, when one calls eq.(l) Hermitian one matrix model, "Hermitian" actually refers to the measure, 

= Yli j=i'^'^i3 (dictated, in its turn, by the norm ||i5$|p = Tr ^$^), but not to the integration contour: 
when one defines partition function as solution of the Virasoro constraints, there is no reason to integrate in 
(1) over the contour = associated with Hermitian matrices. Moreover, since D-module is defined by 
linear equations, the superposition principle is applicable, and a sum of any two solutions is again a solution. 
If matrix integrals like (1) with different N satisfy the same Ward identities like Virasoro constraints (6), then 
linear combinations of integrals with different A'^ still are solutions. This means that the matrix size N can be 
also interpreted as characteristic of integration contour, moreover, this characteristic is not obligatory a positive 
integer, but can be also negative, rational and even complex- valued. 

Differences between possible integral realizations of one and the same partition function are not exhausted 
by freedom in the choice of integration contours, even with extended interpretation of the term "contour". 
Partition function can be represented by absolutely different classes of matrix models. The simplest important 
example is provided by the model (1): it can be represented not only by (1), but also by a Gaussian integral 
[41] from Kontsevich family (3) with tk = ■^trL"'' [42, 6]: 



In this representation analytic continuation to non-integer values of N is even more straightforward. Again, 
linear combinations of solutions with different values of N , are still solutions to the Virasoro constraints (6), 
and all of them should be considered as particular branches of the 1-matrix model partition function. This 
provides additional information about the problem of conjugate variables: N itself can serve as one of them, 
or the coefficients of above-mentioned linear combinations can - changing the former for the latter is a kind of 
Fourier transform in the space of conjugate variables. 

Possibility to represent one and the same partition function by members of two different matrix model 
families, by (1) and by (11), is a typical example of duality between different families of quantum field theories. 
One more duality follows from existence of CFT representations for the same partition function, it relates these 




(11) 
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matrix models to 2d free scalar theory and - through one more duality - to free fermions. These dualities 
illustrate the fact that the nature and the number of integration variables can differ in the most radical way, 
still partition functions coincide: in absolutely different theories one can find families of identical correlators. 
Once again, this is an illustration of the general stringy idea: everything is the same, one can fully describe one 
theory (phase) in terms of another. Matrix models provide a nice framework for testing this principle and can 
help to transform it into constructive and reliable procedure. 



4.1.4 Limiting procedures and asymptotics 

Different realizations of partition functions are adequate or at least convenient for descriptions of different phases 

of the system (branches of partition function), or at least different classes of phases. As usual, phases become 
pronouncely different in particular asymptotics. Systematic approach implies that all kinds of asymptotics 
should be investigated. In particular, in the case of matrix models, all kinds of continuum limits should be 
examined, not only naive or t'Hooft's or double-scaling. Obviously, all kinds of multi-cut solutions and all kinds 
of multi-scaling limits require attention, but this is only the beginning: by no means all possible asymptotics 
are exhausted in this way, and also there is a lot of interesting beyond continuum limits. Unfortunately, not 
too much is known about this variety of limits and nothing like classification of asymptotics is available (this 
was a trivial thing in the case of one variable and becomes a pretty sophisticated issue in the case of infinitely 
many variables). 

There are two different possibilities to perform limiting procedures: at the level of defining equations (linear 
Ward identities or bilinear Hirota-style equations) and at the level of particular realizations. Of course, one and 
the same limit can look very different in different, though equivalent, representations. Because the subject is 
under-investigated, the set of known realizations is poor and many potentially interesting asymptotics cither do 
not attract attention or are difficult to handle. We mention just two old problems, which still remain unresolved. 

In the class of generalized Kontsevich models (3) and (4) one naturally distinguishes [28, 11] between 
Kontsevich phase - the asymptotics of large L - and the characters phase - the asymptotics of small 
L. In these phases the " time- variables" tk = rk + ^trL"*^/" with n = deg W'{z) are respectively close to 
Tfc = fc(„"Lfc) res (VF'(z))^~''^" dz and to infinity. Kontsevich phase corresponds, at least through the relation 
(11), to perturbative phase of the models like (1), while the characters phase is the strong coupling limit of 
these models. In this sense (11) is an example of S-duality, interchanging weak and strong coupling phases. 
Actually, these phases are closely related to the strong and week coupling limits of Yang-Mills theory and thus to 
confinement problem. However, this problem was never exhaustively analyzed even at the level of matrix models 

- despite the existence of Kazakov-Migdal-Kontsevich family (4), which - in variance with realizations like (1) 

- provides an effective tool for studying both asymptotics and interpolate between them (for these models 
Kontsevich phase is the WKB asymptotics, while the characters phase is perturbative limit). A systematic 
analysis, like suggested in [2, 43, 44] for (1) continues to wait for its time for Kontsevich and for more general 
unitary matrix models. 

Another abounded problem would be just a small paragraph in this systematic analysis, still it is quite 
important by itself. This is the problem of "double-scaling" limit [45] of (1), which is one of the simplest 
non-naive large- iV asymptotics of Z(t\N), when N —>■ oo together with certain special adjustment of t-variablcs. 
For investigation of this limit one can make use of various techniques, of which the most important are two: 
taking limit of loop equations (Virasoro constraints) (6) [46] and exploiting the identity (11) and taking limit 
within the family of Kontsevich models [42, 6]. In these ways one can argue that 
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(12) 



where at the l.h.s. t2k+i = 2^+1^^^ ^fe 1 _ q g^j^^j _ -l-trL and the quantity at the r.h.s. (expressed 
through i2k+i = fc-i-i/2 ^^'^~*^~^^^) ^ special name of Kontsevich r-function TK{i) and is an increasingly 
important matrix-model and string-theory special function. However, even reliable derivation of the principally 
important result (12) and exact relation between t, L and t, L is not available, nothing to say about corrections 
to this formula (describing what happens when one approaches the limit) or its generalizations. Of special 
importance among such generalizations are the limits, when Z{t\N) — > rj^it) with n > 1, because of their 
obvious relation to Givental-stylc decomposition formulas [47, 2, 40] and their potential role in applications. It 
goes without saying that various non-perturbative asymptotics of tk should be investigated: this is the next 
natural task after such program is put on track for the basic special function Z{t\N). 
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4.2 Integrability 

Integrability is believed [1, 6] to be the pertinent property of quantum partition functions, reflecting the fact 
that they are results of integration, which erases most of initial information (roughly speaking, almost all forms 
arc exact), leaving only cohomological variables. Therefore integrability is intimately related [48] to topological- 
cohomological-holographic-stringy theories and all these kinds of ideas should be considered together. However, 
despite many efforts, spent on investigation of particular models, we are still far from formulating clear concepts. 
Partly this is because examples are dictated more by applications than by the internal logic of the theory, and 
they provide somewhat chaotic flow of information. There are already several things that we seem to know, but 
their complete list and relative significance remain obscure. 

4.2.1 Hirota equations 

The first in the list are bilinear Hirota-style equations, which establish contact with Lie algebra theory [27]. 
Bilinearity is usually related to the properties of comultiplication and to the basic relation A{g) = g ^ g for 
group elements [49, 50]. Connes-Kreimer theory [50]- [52] implies, that quantum partition functions respect 
these properties, even if they are defined through Feynman diagrams, without explicit reference to functional 
integrals, and identify the underlying algebraic structure with (generalized) renormalization group. As already 
mentioned, deeper understanding of relation between integrability and renormalization groups remains among 
the main open problems of the theory. 

4.2.2 Moduli space of solutions 

The second big problem is midcrstanding, classification and control over the freedom in solving Hirota equations. 
We mentioned already, that these bilinear equations restrict the dependence on coupling constants tk, and the 
freedom remains in the conjugate/dual variables, like zero-modes, or boundary conditions, or choices of vacua, 
or holographic data, or whatever else name and analogy one prefers to use in connection with one's favorite 
application. The problem is to find not just a name, but adequate language to speak about this data. Suggestive 
example is provided by original Hirota equations [53], associated with KP hierarchy: there the freedom (or at 
least a part of it) can be interpreted in terms of Riemann surfaces. Different KP r-functions, at least from the 
class of the finite-zone solutions, differ by the choice of Krichever data E [54] : a complex curve, a point on it 
and holomorphic coordinates in the vicinity of a point. Given this data, the KP r-function is fixed to be 

T(t|S)=e*«*0(^Bfctfclry), (13) 

k 

where is the period matrix of the Riemann surface, Bk = §g^k and ^k{0 = (^~('^+^) + 0{1)) d( are 
meromorphic 1-forms with vanishing ^-periods. If ansatz (13) is substituted into Hirota equations, they become 
equivalent to Shottky condition for the period matrix Tij [55], i.e. require 9{z\Tij) to be Riemannian thcta 
function, not just arbitrary Abelian one. Rather straightforwardly, the simplest - rational and solitonic - KP 
T-functions can be treated as particular cases of (13) with Riemann surface of genus zero. All remaining, 
non-finite-zonc KP T-functions are more sophisticated and can be considered as associated with infinite-genus 
Riemann surfaces, though exact meaning is not yet ascribed to these words. Clearly, there are many different 
kinds of "infinite-genus" T-functions. Even the principal example of Kontsevich t- function Txit), which belongs 
to this infinite-genus family, is not well understood and investigated. 

Anyhow, example (13) remains the main reference point in the theory of integrable systems. It underlies 
the belief that the right language for description of dual/conjugate variables to times (coupling constants) 
should be that of Hodge theory, different T-functions should be parameterized by data, encoded in spectral 
surfaces, not obligatory 2-dimensional. Thus integrability is believed to be a unification of group theory and 
algebraic geometry, with certain combinatorial flavor, already coming from the studies of adjacent and clearly 
related subjects [34, 50, 51]. At the same time, even despite these ambitious perspectives, too much remains 
obscure in the general structure of integrability theory, see [56]- [71] for various important constituents which it 
should finally uniiy. All this makes natural the fastly increasing attention to integrability concepts, but decisive 
conceptual breakthrough is still to come. 

4.2.3 Seiberg-Witten theory, Whitham integrability and WDVV equations 

The third problem is the search for an adequate language, which can help to merge Hodge theory with Lie algebra 
structures. There are various approaches. The most obvious is to study cohomological (topological) field and 
string theories [72] and find traces of integrability there. This program is very successful in dealing with various 
examples, but general concepts appear very hard to extract from it. Among such concepts definitely are the 
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WDVV equations [73, 74] and Batalin-Vilkovisky formalism [75, 76]. Alternative approach, from integrability 
side usually starts from (13), directly or implicitly, and notes that theta-function in that formula is oscillating, 
therefore one can distinguish between fast (oscillating) and slow variables. Then, as usual, one can take an 
average over fast variables and consider an effective theory induced on the space of slow variables. Since, 
according to (13), the fast variables are exactly the times tk, the slow variables of the emerging effective theory 
should be exactly the dual variables, that we are interested in. This idea, known as Whitham theory, has 
obvious contacts with those of renormalization and renormalization group, and not-surprisingly appears to have 
immediate applications. While absolutely un-developed in general, it has a solid very well formulated chapter: 
the Krichever-Hitchin-Seiberg-Witten theory, describing explicit construction of peculiar special functions — 
Seiberg-Witten prepotentials - from the analogue of Krichever data, associated with peculiar Seiberg- 
Witten families of spectral surfaces (the best known examples are Calabi-Yau spaces and special families of 
Riemann surfaces). The prepotential J^{Sk) is defined [77]-[79] by two equations: 



where flswiz) is a linear combination of f2fc(z) in (13) and holomorphic 1-differentials oj{z) (with non-vanishing 
^-periods), possessing a special- geom,etry property that SQsw is a symplectic form on the "sheaf of spectral 
surfaces over the moduli space. In practice this condition states that 6Qsw is holomorphic (has no poles) if the 
variation is taken along the moduli space of Seiberg-Witten curves. Canonical system of contours {Ak,Bk} in 
(14) consists of all non-contractable cycles (including those around resolved singularities of Ofc) and their duals, 
see [79] for details. The variables Si define a sort of flat coordinates on the moduli space. 

It appears that the so defined prepotential always satisfies the (appropriately generalized) WDW equations 

[74], namely for matrices (^ij = Qs-fs^ds ^^'^^ out of prepotential's third derivatives. 



for any triple (z, j, k). Moreover, this system of equations for infinitely-large matrices often survives reduction 
to finite-dimensional matrices when only S'-variables, associated with ordinary non-contractable contours, are 
taken into account. This reasons why such reductions exist are far from obvious, and so is the entire theory 
of WDVV equations. Technically, it relies upon residue formulas for moduli derivatives of period matrices, 
but conceptually it involves reduction theory of non-trivial algebra of forms and is still very unsatisfactory. 
Additional puzzle is that reduced WDW equations are violated, at least when applied naively, in the case of 
elliptic Calogero system (see the third paper in ref.[74]). Other elliptic examples [80, 71] were not yet analyzed 
in the context of WDVV equations. 

At the same time the theory of WDVV equations can appear even more important than it seems. Today this 
set of equations is the only candidate for a definition of a prepotential in internal terms, which does not refer to 
explicit construction procedure. Thus WDVV equations can probably be promoted to the status of definition of 
the prepotential and acquire the status, similar to the one that Hirota equations have for r-functions. (Referring 
to its possible origin as a Whitham average of (13) - or, more, accurately, of the underlying Hirota equations,- 
prepotential is often called Whitham or quasiclassical r-function, and the above-mentioned problem is to find 
the miernaZ definition of quasiclassical integrability - a very important theoretical challenge.) The main problem 
with WDVV equations from this point of view is that they describe only " spherical prepotentials" : from the 
study of topological theories and related algebraic constructions [81] we know that there is a whole hierarchy of 
prepotentials, associated with the genus expansion of Yang-Mills and matrix models partition functions, and 
Seiberg-Witten prepotential is just the zeroth (genus-zero) term of this sequence. Appropriate generalization of 
WDVV equations for entire sequence is not yet found, for first steps in the cases of genus one and two see [82], 
and for relation to Batalin-Vilkovisky theory see [83] (to avoid possible confusion, these papers consider only 
WDVV equations with additional requirement of " unit metric" , which is natural in applications to quantum 
cohomologies [84] and naive topological models - to " geometric prepotentials" ,- but is not quite straightforward 
to relax). 

4.3 CFT representations, Wick theorem and decomposition formulas 

The theory of 2d conformal systems [85], once at the front-line of attention in string theory, was suddenly 
abandoned before an exhaustive theory was formulated, and many problems, including effective description of 

entire variety of conformal models and equivalencies (dualities) betwcicn them, are left unresolved. Partly, this 
happened because attempts to glue the entire subject together - in the framework of landscape theory [35] - 




(14) 
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naturally caused a desire to interpolate between different conformal models, and this unavoidably embeds [32] 
the world of 2d conformal theories into that of 2d integrable systems - still a very badly understood subject, 
intimately related to the theory of quantum Kac-Moody algebras [86, 87, 88], which attracts much less attention 
than it deserves. 

As already mentioned, there is no doubt that every class of matrix models have CFT representations, i.e. 

partition function of a given family of matrix models can be represented as a correlator of some operator in 2c? 
conformal theory. The questions are: what is the way to build this operator - so far it is matter of art rather 
than a systematic procedure;- what kind of information from the much bigger phase space of 2d free fields is lost 
(projected out) in transition to the matrix models (in particular, even models which are not dual themselves 
can still produce representations of one and the same matrix model partition function); how the data of matrix 
model is mapped into the data of conformal model etc. Today we do not know anything similar to answers to 
any of these questions. 

The crucial property of free field (and thus of 2d conformal) theories is straightforward realization of the 
Wick theorem: all correlators are decomposed into multi-linear combinations of pair correlators (see [89] for a 
rather fresh discussion of the issue). In principle, for topological and integrable theories naive Wick theorem is 
not applicable: even after all possible simplifications there are contact terms (and prepotential J-{S) is almost 
never quadratic). In interacting theory transformation to angle-action variables which would explain to what 
kind of correlators the naive Wick theorem is applicable, is highly non-linear and probably not very useful. 
A very important question is what substitutes the naive Wick theorem for generic partition functions. As we 
know, generic partition functions are not at all in general position in the space of all functions, they satisfy 
many constraints, like linear Ward identities, bilinear Hirota equations, highly non-linear WDVV equations. 
Thus we can expect that above question about Wick theorem can make sense. Of course, one can try to 
interpret as avatars of the Wick theorem all of above-mentioned relations on partition functions. Remarkably, 
there is also a much more profound candidate: decomposition formulas. It appears that partition functions 
can often be represented as poly-linear combinations of simpler building blocks, like Gaussian matrix models 
and Kontsevich r- functions tk- So far this fact, which we call Givental-style decomposition, was explicitly 
formulated in restricted number of cases [47, 2, 40], but we believe that this is just the tip of the iceberg, 
and further work will confirm the universality and important role of decomposition phenomenon: a probable 
substitute of Wick theorem for generic partition functions. 

5 Towards exhaustive theory of 1-matrix model: the Dijkgraaf-Vafa 
theory 

Dijkgraaf-Vafa theory [90]-[104] deserves special attention, because for the first time an application was fomid, 
which requires understanding of a whole variety of phases of a single matrix model, namely the model (1). 
Despite it is still not the entire variety - only peculiar phases are considered interesting, where partition function 
is consistent with the genus expansion, - it is still a very big step forward, and it stimulated fast progress in 
matrix model theory. In this section we briefly characterize some directions of this progress As everywhere in 
these notes, we concentrate on pure theoretical issues and ignore the - sometime very interesting - results in 
applications to Yang-Mills, quantum gravity and model building. 

5.1 Genus expansion 

Genus expansion, i.e. the t'Hooft's 1/iV expansion [105], attracts constant attention since the early studies of 
matrix models [15]. While it has clear meanings both in diagram technique (genus expansion for fat graphs) and 
in WKB approach to the integrals like (1), specification of partition functions possessing genus expansion among 
generic solutions of Ward identities [17] is not straightforward and looks somewhat artificial. As explained in 
[2, 44], it involves several steps. 

First of all, by rescaling of all time- variables tk ^tk one introduces the new parameter g. It appears in 
Virasoro constraints (6) as a coefficient in front of the double-derivative terms. Second, solution of these 
Ward identities is looked for among rather special functions, such that 

5^ log Z{tlg\N) = J2 g^'F^"^ {t\9N) (16) 

p>0 

is a series in non-negative powers of g^, provided the coefficients F^P\t\gN) of this expansion depend on the 
t'Hooft's coupling constant S = gN . Though consistent with (6), this is actually a tricky requirement. It 
implies that Z itself is a series in all integer, positive and negative, powers of g^, just coefficients in front of 
negative powers are strongly correlated. Moreover, there is no N in Ward identities (6). N can be introduced 
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as S* = gN = d{g^ logZ)/dto, but this requirement breaks superposition principle: once it is imposed, the sum 
of two solutions is not a solution (this is nearly obvious: a sum of two exponentials is not an exponential). This 
means that genus expansion can be at best a property of the elements of a linear basis in the space of solutions, 
while generic element of this space does not respect it. In other words, genus expansion can be a property of 
particular integrals, like (1), but there is no way to impose this requirement on linear combinations of integrals 
with different N, which will still be solutions of the Virasoro constraints. 

There are at least two ways to deal with this problem (the obvious possibility to ignore it, as is usually 
done, is not counted). First, one can look for the ways to define special basises - the problem is to define them 
in invariant way, by specifying their properties, not by explicit construction, like matrix integrals. Second, one 
can try to formulate genus expansion without explicit reference to N - this road leads to introduction of check 
operators [44]. 



5.2 Gaussian and non-Gaussian partition functions 

Requirement of genus expansion is not enough to specify the phase/branch of partition function completely, it 
regulates only dependence on the common scaling factor of all couplings tk and says nothing about their ratios. 
The next step [2] is to select combinations of ^-variables which can appear in denominators. A way to do this 
is to make a shift, tk ^ tk — Tk, and consider W{(j)) = J2k'^k4''' as a bare action and as perturbations. 
In other words, the branch 

= VoK^ I./^^^P [-TrWm+Y.*kTr (17) 

of partition function is defined as a formal series in non-negative powers of tk, while W{z), i.e. T^'s, are allowed 
to appear in denominator. Once W{z) is introduced, one immediately distinguishes between the Gaussian (with 
W{(j)) = ^^(f^) and non-Gaussian phases. 

Next, it turns out that while the variable N is naturally introduced in the Gaussian phase, in non-Gaussian 
case the situation is more sophisticated: there are naturally n = deg W'{z) parameters like N: the phase is still 
split into a n-parametric family of phases. Additional variables can be interpreted as numbers of eigenvalues 
concentrated near the n extrema (maxima and minima) of the bare potential W{4)), a more reliable interpretation 
is as integration constants in solution of shifted (by W{z)) Virasoro constraints. See [2] for further details. 

Gaussian partition function can be represented as 

zi'\t\S\g) = exp (1 (-5To - ^5^)) exp {Y.g'^-'T^i\t\S)^ (18) 

and information about the t-dependence of prepotentials TQ\t\S) can be represented in terms of multi- 
densities, the generating functions of m-point correlators, which - for given p and m - are ordinary poly- 
differentials on the Riemann sphere. 

Multi-densities can be defined in different ways (for different families of correlators) [2], one of the ways 
consistent with genus expansion and with the free-field representation (6) of Virasoro constraints makes use of 
the operator V{z) = 2(i<^_(z) = '}2t=o — — ' 



V(^i) . . . V(^„) (g^ log Zw{t\S\g))] = p^-\^\zi, z^\S\g) = ^sV^''"^^!, • • • ,^m|^) (19) 
The first Gaussian multi-densities are: 



p>0 
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We assumed here that M = 1, then yQ = z"^ — Av and v = S = gN . 

Non-Gaussian are somewhat more sophisticated, they are actually made from Riemann theta-functions for 
a pecuhar Seiberg-Witten family of hyperelliptic complex curves. For example, the (p, m) = (0, 2)-density 
is the 2-point correlator on such surfaces: 

Pw'^\zi,Z2) = d^,d^^ \ogE{zi,Z2), (21) 

where E{zi,Z2) ~ "'elzl-'z^)'^ is the prime form [5, 3]. 

5.3 Decomposition formulas 

Rewriting (17) in terms of the eigenvalues and then separating the eigenvalues into sets, associated with different 
extrema of W{(j)), one can deduce [93] a decomposition formula [2] for non-Gaussian Zw- 



wo\{u{Ny) 



n+1 



n^^^ e-^^^("^)Vol(t/(jV.)) fr 2N,N,^ fy. TT.(M.).,W|^. .,2^ 

with ai denoting the n roots of 

W'{z) = J2 kTkZ^-^ = {n + l)T„+i JJ(^ - a^), 
Oiij = ai — aj, the frequencies Mj = W"{ai) = (n + l)T„+i Ilj^i operators 

a,=expf2f:( -(^ + ^-1)' ' ' 



= exp 



fe>3 "^k 



As explained in the previous section, this Givcntal-stylc formula is an example of general phenomenon, gener- 
alizing the Wick theorem from free fields and conformal theories to matrix models and non-trivial integrable 
systems. 

Eq.(22) requires (and deserves) deep investigation. Today almost nothing is known about it. It is unclear 
how (22) follows directly from Virasoro constraints, without explicit use of matrix integrals. It is unclear what 
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is exact implication of the relation (12) for Gaussian functions, standing at the r.h.s. of (22). It is clear that 
the combination of (12) and (22) should represent Zw as some operator acting on r^, but exact formula is 
unavailable (Greg Moore in [39] and especially Ivan Rostov in [40] came very close to the answer, but it still 
escapes). 

Eq.(22) can be also considered as a definition of a linear basis in the space of solutions of Virasoro constraints, 
which is formed by the functions, respecting the genus expansion. Arbitrary solution can be represented as a 
linear combination, 

Zw{t)= [ fi{S)Zw{t\S), (23) 
Js 

with arbitrary measure i^{S). 



5.4 Check-operators and t-evolution operator 

According to their definition (19), multi-densities do not depend on the time variables tk, only on Tg, . . . ,Tn+i 
and S. Actually, dependence on T„ and T„_|_i is fixed by (6), and multi-densities depend on 2n variables: n 
of them are T's and n are S"s. Therefore one can ask for a procedure which defines multi-densities directly in 
terms of these variables, without an intermediate introduction of - infinitely many - auxiliary (from the point 
of view of multi-densities theory) variables tk. This problem is solved, at least conceptually, in [2, 44] in terms 
of check-operators. 

Multi-densities are not arbitrary functions of their 2n variables, the Ward identities (6) express all of them 
in terms of a single function of n -I- 1 variables. To, ... , Tn-i and g [2], which we call bare prepotential. The 
choice of bare prepotential is actually the choice of particular branch of partition function. Note, that this 
means that there are not just infinitely many branches, they form a continuous variety and this is a property 
of every reasonable partition function and of anything obtained by application of evolution operators (for 
which the functional integral is a particular realization), sec [34] for explanation, how continuous phase structure 
emerges from discrete bifurcations (phase transitions), and for discussion of the adequate formalism. This set 
of questions is very important, in particular it can help to understand the relevant topology and may be even 
the metric structure on the space of phases, what is the principal problem in landscape theory. Development of 
such formalism for matrix models is a task for the future, today we know how to proceed within a given phase, 
when bare prepotential is somehow specified and the question of relative importance of different choices is not 
addressed. 



5.4.1 Independent variables and bare prepotential 

To understand what it means to specify the bare prepotential, we consider immediate corollary of (6) with 
tk — > i(ffe — Tfc) for partition function at vanishing times tk- It is given by [44]: 



Z{T\9)\t=o = J dkz{k\m, r?n|/)e*^'"-''"^ (24) 
with an arbitrary function z of n arguments {k,r]2, - ■ ■ , r]n) and g^. Here the Z/_i-invariant variables are used, 
«^ = ^logT„+i, r,k=T~^ {T^ + ...), x = To + ...r^r,n+i (25) 

In particular, 



„ _ {rj.2 2(n + l) \ 



V3 - y-in J-n-l-Ln-Ln+1 + -'-n-2-Ln+l J ' 

4(n+l)^ 8(n+l)2 8(n + 1) 



3n 



( k , fc(fc-2)! y^ , (n+l)'(n-0! k-i-ir^i \ j.-^ 

The variable x is obtained from rjn+i by normalization and it is the only variable which contains Tq. At the 

same time only To appears in double-derivative item of Lq, thus in (24) the io-constraint links the x- and 
w-dependencies. Of course, generic (24) does not possess genus expansion, i.e. does not guarantee that the 
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bare prepotential F{T\g) = g'^\ogZ(T\g), is expanded in non-negative powers oi g^: as explained above, this 
requirement is somewhat artificial from the perspective of D-module theory. Still, a simple ansatz makes things 
consistent (it is not quite clear if this ansatz is absolutely necessary for genus expansion to exist): if 

dF 

S= — = const (26) 

i.e. is independent of Tq, . . . , T„+i and g, then, 

z{k\v2, Vn\9^) = H{r]2, Vn\S\g^) S{k - S) (27) 

and ^ 

F = gHogZ = Sx + ^ logT„+i + g^ log H{r,2, Vn\S\g^) (28) 

where H is an arbitrary function of n + 1 variables. Genus expansion occurs, if we request g"^ log H is expanded 
in non-negative powers of g^. 

Introduction of the other S^- variables can be considered as ingenious version of Fourier-Radon transform 
with the help of Dijkgraff-Vafa partition functions (22), which converts H{ri2, . . . ,r]n) into an arbitrary measure 
H{S) of the 5-variables in (23) [2]. 

5.4.2 Check-operator multi-densities 

Other L/j-constraints with fc > express ^-dependencies of Z[t\T) through T-dependencies. In particular, 
they allow to represent multi-densities (both Gaussian and non-Gaussian) by action of T-dependent opera- 
tors on partition function (24). Such operators, involving only derivatives with respect to T- variables, arc 
named check- operators in [44], to distinguish them from hat-operators, acting on t- variables. The task is to 
express multi-densities, defined by application of hat-operators to the full prepotential F{t\T\g), through check- 
operators applied to the bare prepotential F{T) = F{t = 0|T|g). These check operators will predictably be 
more sophisticated than hat-operators ^{z), but instead they can be applied to the independent (free) function 
F{T). 

Remarkably, the problem of building check-operator multi-densities appears essentially equivalent to the 
problem of Gaussian multi- densities [44]. Though this fact is not fully proved - and even adequately formulated 
- yet, it is not too surprising: the needed check-operators are universal, in certain sense they do not depend on 
the phase, thus they should be restorable from information, available in any phase, including the Gaussian one 
(once again we encounter the main idea of string theory). Anyhow, modulo some details concerning ordering 
prescriptions, the first check-operators can be just read from formulas (20) for Gaussian multi-densities [44]: 
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Here 

f{z\g) = W'{zf+4g^Rw{z) and (z) = ^ (a + 6 + 2)i;+b+2z" 

a, 6=0 

The check-multidensities are defined to satisfy 

^^"'(zi, . . .,z^\g)Zw{T\g) = A'|^"'(zi, . . . , z^\g)Zw{T\g) = 
= V{zi)...^{zm)Zwmg)\\^^^ 
for full correlators Kly"^\zi, .. .,Zm\9), which are derivatives of Z{t), not of its logarithm, while for connected 

Pw"'\zi, ■ ■ ■ , Zm\g) 

{Zw{T\g)r^ p'itHzu Zm\9) Zw{T\g) + • • • , Zm\d) = 

= Viz,)...Viz„,){gHogZw{T\g)) 

Relation between the fuU and connected correlators is provided by straightforward - though heavily looking - 
combinatorial formula [44] 
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(31) 



which serves as generalization of the Wick formula from the case when the only connected correlators are 2-point 

(note that this simplification docs not occur even in Gaussian phases of matrix models beyond naive continuum 
limit). A similar expression relates the check-operators K and p, in particular, 



K^4"\^-,9) = Y.9''-'pT{^-,9) 

p=0 



(32) 



Note also, that in variance with (20), expressions in (29) explicitly contain g"^. 

Actually, once the check-operators multi-densities are introduced, one can forget about the genus-expansion 
constraint: operators can be applied to any bare partition function. Instead of being a constraint on partition 
functions, genus expansion requirement dictates becomes a selection criterion for a basis in check-operators space. 
The next step is to build up a theory of the t-evolution operator U{t), which generates the f-dependence of 
partition function: 

Z{t\T) = IJ{t)Z{0\T). (33) 
Its existence is almost obvious, see [2, 44], but properties and explicit realizations remain to be found. 

5.5 Seiberg-Witten theory 

The Dijkgraaf-Vafa theory is intimately related to Seiberg-Witten theory. 



15 



5.5.1 Genus-zero level 



For detailed discussion of this relation in the particular case of Dijkgraaf-Vafa partition functions per se, i.e. 
with requirements that genus expansion exists and only genus-zero contributions (spherical prepotentials) are 
considered, see [96]-[100]. Here arc the main points of this analysis. The smallest moduli space of Sciberg- 
Witten structure unifies all phases with all the bare potentials W{z) of a fixed degree n. It has (complex) 
dimension 2n, T and S'-variables form the set of 2n flat moduli of Seiberg-Witten structure, defined with the 
help of risw = yJW'{zf - ^jw{z)dz, where fw{z) = -Rw{z)F'-"\T) is a polynomial of degree n — 1, made 
from arbitrary function F^^\T) of n variables To,...,T„_i (and dependence on Tn and Tn+i prescribed by 
L-i and Lo-constraints). The space of all functions of n variables can - in the framework of Fourier- Radon 
transforms - be parameterized by arbitrary functions of n other parameters, and S can play the role of such 
parameters no better, no worse than any other set of variables. What distinguishes S is that they are flat 
Seiberg-Witten moduli, i.e. the symplectic form on the sheaf of Seiberg-Witten curves is 

Silsw ^i^tiz) A6Si + nk{z) A6Tk (34) 

and 

5^ = (jli^flsw, Tk = lesaoz'^ilsw (35) 
The CIV-DV prepotential [90, 91] is defined by 



d^civ f ^ dTciv 1 _fep, 

= * VLsw, — = 7reSoo-2; ^sw (36) 



dS Is ' dTk k 

it can be explicitly calculated term by term and possesses pronounced group-theoretical structure [98, 2], which 
still awaits its adequate interpretation and explanation. No reasonable generic formulas for this prepotential are 
yet found, despite it is a genus-zero quantity, and despite it can be deduced in at least two dual ways: from (36) 
and from (22). Moreover, even direct relation between these equivalent representations is not yet established. 
Also important is that there are various differently looking forms for the same items in J-ciVi their equivalence 
follows from peculiar identities [99] , which are simple to prove, but not so simple to understand. 



5.5.2 WDVV equations 

In [96, 99, 100] the proof is given of the WDVV equations (15) for the CIV-DV prepotential. This is an additional 
check that 2n moduli T and S form complete set, at least at the level of genus-zero prepotentials. In [97] a less 
obvious observation is reported: sometime this set can be reduced further, for example, by elimination of one 
of the 4 moduli in the case of n = 2, however, the meaning of this observation remains obscure. 

The analogues of the CIV-DV prepotential for higher genera can also be defined because of its relation to 
Ward identities (6) - this makes the situation different from generic Scibcrg-Wittcn theory, where the problem 
of lifting from genus zero to all genera remains very important, but unsolved. This makes it possible to ask, 
if higher genus prepotentials satisfy the (properly generalized) WDW equations, for example the ones implied 
by studies of topological models [82]. This problem was not seriously attacked yet. 



5.5.3 Towards quantum Seiberg-Witten theory of check multi-densities 

The next natural question in the theory of check operators concerns definition of S- variables. Since in particular 
phases they are periods of P^y)^\ it is natural to treat 5- variables as operators, defined as periods of p^^l^^ 
Moreover, under certain conditions the A-periods of higher-genus one-point densities are vanishing, and one can 
even consider the periods of p^ l^^ or 

Si = £-J^ K^-^'\z). (37) 

This is a very perspective direction of research. Among other, obvious and not-so-obvious, questions is devel- 
opement of the full-scale version of Seiberg-Witten theory for check-operator quantities, i.e. a quantization of 
Seiberg-Witten structure. Introduction of check operators with the goal to deal with arbitrary bare prepoten- 
tials in a uniform way is a nice illustration of string "third-quantization" procedure, which implies that in order 
to study families of models one actually needs to quantize a representative of the family. An example of the 
simplest relation in this quantum Seiberg-Witten theory is provided by the quantum version of (14) [44]: 



(38) 
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probably, with no corrections at the r.h.s. 

The check-operator A-and B-periods in (38) are equal respectively to Si = S^^^ + g'^S^^^ + 0{g'^) and 
nf^ +0(^2), where [98] 



(39) 



2^(1) (29) W"'{z)dz _ d^Wd^W - {d^Wf 
9 - res„, g^,^^-|2 - 8(92 W^)3 

and ai, i = 1, . . . , n are the roots of W'{z), operators Ri = Rw{oii), 



M, = W'iai) = {n+ l)T„+i [] ay, 

Uij = ai — aj. Note that in the leading order the operator B-periods and genus- 1 contributions to A-periods 

do not containt T-dcrivativcs, thus non-trivial arc only commutators, involving Sf'^. These commutators can 
be easily deduced from two general formulas from [44], 



[Rw{x),W{z)] = 



^M^^, [^.(.),^.(.)] =(a.-a.)^kM^^, 



(40) 



together with 



k-l 



Mi ' 



from [99], which implies straightforwardly: 

[Rw{x), a,] = ^ (a + b + 2)Ta+h+2x'' = -— !, _ ^ 



a, 6=0 



and, since W'{ai) = W'{aj) = 0, 



= [Rw{ai),aj] = — 



(41) 



(42) 



Eqs.(40) need to be supplemented by (42), because in commutators of periods (39) the arguments at also depend 
on Tfc and their derivatives should be taken into account. Thus instead of (40) we need: 



[Ri,Mj] = [Rw{ai),W"{aj)] = 8: 



W'{ai) - W'{z) 



ai - z 



+ W"'{aj) [Ri,aj]=- 



2Mi 



a: 



3 ' 



(43) 



and 



[Ri,Rj] = [Rw{ai),Rw{aj)] = {dz - 9x) 



Rw{x) — Rw{z) 



X — z 



+ 



+ [Rw{ai),aj] Rwioij) - [Rw{aj),ai] Rw{ai) = [Ri - Rj) 



Now we have everything to check, that 



^(0)^^(0) 



9 1 



9 



MiMj 



[Ri, Rj] -J^[Ri^ Mj] Rj + — [Rj,Mi\ i?, ) = 



Mi 



Direct generalization of (43) states, that 



\Ri,d"'W(ai)\ = -m\ y — —7 j-^ 



fc=0 



in particular. 



so that 



Ri, /Q2M/\3 ^^i> 



{d'^W) 



MiMjal 



24 



(44) 
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is symmetric under the permutation i j. Since also g"^ s'l^\ g"^ S^'^^ 
obtain: 



and 



g's^\tif' 



0, we finally 



fT(0) Yi(0) 
j 



n 



[W{ai),W{aj)]=0, 



Sl^+g^s^\ 



n 



(0) 



r „2 



W'{ai)-W'{aj) 



W"{ai) 



g^6i, 



(45) 



Note that in the last formula there is no need to vary the argument aj, because this variation gets multiplied 
by W'iaj) = 0. When j ~ i one needs to apply I'Hopital rule to resolve the 0/0 ambiguity, and it provides the 
non- vanishing answer. This ends the proof of (38). It is an open question, if there are higher order corrections 
~ 0{g^) to the r.h.s. of that formula. 

Given (38), of special interest is this kind of quantization of the WDVV equations. At least in principle, 
such approach opens a possibility to unify WDVV equations for different genera, the spherical equations from 
[73, 74] with higher genus equations from [82, 83], and thus provide a new deep connection with geometry of 
the moduli space of punctured Riemann surfaces. 



6 Neoclassical period: back to concrete results and back to physics? 

The story about the challenges of matrix models is nearly infinite and can be continued further and further. 
We, however, draw a line here. The only brief comment that deserves being made at the end of these notes, is 
that the transcendental period, when most attention is concentrated on non-obvious structures underlying the 
subject and thus on high abstractions, is never the final stage of dcvelopement of a theory. Sooner or later, 
deep insights from transcendental studies open a possibility to solve practical problems and, in our case, derive 
concrete formulas with the left hand side and the right hand side. Remarkably, some of such results are already 
starting to emerge. Despite they are rather week, it is principally important that they appear, it means that 
the theory is on the right track and new knowledge continues to have many non-trivial intersections with the 
old - what is a necessary feature of a healthy research project. 

We mention here just three results of this kind, all concern the theory of CIV-DV prepotentials and all 
establish non-trivial relations with free fields on Riemann surfaces: the previous-level theory, from which one 
starts from when departing into the world of integrability and matrix models. The first result [104] identifies 1- 
loop correction to the CIV-DV prepotential with certain combination of determinants in some yet unidentified 
- conformal theory. The second result [40] provides additional evidence in a form of a CFT representation 
(not fully polished yet) of the all-genus matrix- model partition function Zw As expected, this representation 
actually constructs the partition function by a Givental-stylc procedure from 2n [n ~ dcg W'{z)) Kontsevich 
fau- functions tk, obtained by the application of star operators [39] at 2n points of the bare spectral curve 
(Riemann sphere), which become ramification points of the hyperelliptic spectral curve, associated with Zw 
The third result [43] formulates a kind of a perturbation theory for correlation functions of 2n operators inserted 
at these points and represents the partition function as a combinatorial " quantum field theory" , essentially of 
Chern-Simons or, better, Batalin-Vilkovisky type. We consider this small set of preliminary results as a clear 
manifestation that all the ideas, discussed during the transcendental period can indeed be brought together and 
the head can catch the tail: the seemingly transcendental ideas will finally form a dense network of knowledge, 
not a shaky road leading into nowhere... 
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